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For a charged particle in a homogeneous magnetic field, we construct stationary squeezed states 
which are eigenfunctions of the Hamiltonian and the non-Hermitian operator X<s> = X cos &+Y sin $, 
X and Y being the coordinates of the Larmor circle center and $ is a complex parameter. In the 
family of the squeezed states, the quantum uncertainty in the Larmor circle position is minimal. The 
wave functions of the squeezed states in the coordinate representation are found and their properties 
are discussed. Besides, for arbitrary gauge of the vector potential we derive the symmetry operators 
of translations and rotations. 

x) : 

The problem of a charged quantum particle moving on a plane subject to a homogeneous magnetic field is met in 
various physical contexts, and it has been extensively studied in the literature and presented in text-books (lj. A 
specific feature of the problem is that the energy spectrum given by discrete Landau levels is multiple degenerate: 
the number of linearly independent eigenstates belonging to the iV-th Landau level is proportional to the area of the 
plane accessible to the particle. The degeneracy is related to the translational invariance: As a classical Larmor circle 
' can be put anywhere in the plane, a suitably defined operator of translation T a = e^ a ' p , P being the generator of 
£ — , translations, commutes with the Hamiltonian and upon acting on an energy eigenfunction tpN(f) produces another 
eigenstate, shifted in space: \T a ^ n {v)\ 2 = n (r ■ + a)\ 2 [^|j3|. Non-collinear translations do not commute in a magnetic 
field, and [P x , P y ] ^ 0. The existence of two non-commuting Hermitian operators, P x and P y , each of them commuting 
with the Hamiltonian leads jlj to the degeneracy. 

The stationary wave functions corresponding to a degenerate Landau level may be chosen to be eigenfunctions 
of either P x or P y (but not both simultaneously). The eigenvalues of P x y are real, and the translation operator 

e iap * (or e lap y ) applied to the corresponding eigenfunction gives only an overall phase factor. The modulus remains 
unchanged by the translations, so that the eigenstates must be infinitely extended in the x— (or y— ) direction. The 
wave functions of an electron in a magnetic field first found by Landau [0J^] give an example: factorized as e tpx x{v)i 
they are eigenfunctions of P x and are infinitely extended in the x— direction (strip- like states). 
As discussed later, in the relation 
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P = -BxR, (1) 



R = (X, Y) has the meaning of the operator corresponding to the classical coordinate of the Larmor circle center 
(the guiding center); X = -^P y and Y = —-^P x . The variable R has a simple classical interpretation, and for this 
reason, it will be used below rather than P. 

Instead of X or Y (P x or P y ), one may choose their Hermitian linear combination Xj> = Acos$ + Fsin<I> with a 
real $. The corresponding eigenstates are "strips" the orientation of which depends on the angle A different class 
of states can be obtained if the wave function is chosen to be an eigenfunction of the non-Hermitian operator A"$ with 
a complex "angle" $ = $i + i$2- By virtue of the relation in Eq. (jl|) , the eigenfunctions of X<$, are also eigenfunctions 
of P®+^ ■ In the case of a general complex <I>, eigenvalues of the non-Hermitian operator P$>+^ are complex numbers, 
and the above argument concerning an infinite extension of the state is not applicable; the eigenfunctions turn out to 
be localized {i.e. the wave function vanishes at infinity). 

In the terminology of quantum optics (for a review see || and references therein) these states belong to the class of 
squeezed states, generalization of the coherent states. In optics the squeezed state is defined as an eigenfunction of a 
non-Hermitian operator, x — iXp, built of two non-commuting variables, the coordinate and momentum of a harmonic 
oscillator (A being the squeeze parameter). A distinctive feature of squeezed states is that the quantum uncertainty 
in the non-commuting variables, is as minimal as allowed by the uncertainty relation (minimum uncertainty states). 
The purpose of the paper to analyze properties of the squeezed states, eigenfunctions of 

Solutions to the Schrodinger equation for a charge in a magnetic field corresponding to non-spreading wave packets 
with a classical dynamics - the coherent states in the modern terminology - were first built by Darwin as early as in 
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1928 0]. More recently, the coherent states in the magnetic field problem have been extensively studied by Malkin 
and Man'ko jjj, and Feldman and Kahn ]§| (see also Ref. J§-|ll[]). In the coherent states, the quantum uncertainties 
in the X and Y coordinates of the Larmor center are equal. Various generalizations to the squeezed states have 
been done by Dodonov et al. ||l2] , ^3| and Aragone In the general squeezed state, the uncertainty in one of the 
coordinates is reduced at the expense of the other one so that their product remains intact. 

In the present paper, we consider stationary states, building the squeezed states from the energy eigenfunction 
belonging to a given Landau level. Being stationary, these states are of different kind than the moving squeezed wave 
packets of Ref. @@. 

The paper is organized as follows. In Section |, we review some general features of the quantum problem of a 
charge in a magnetic filed. In Section O, we define the squeezed states and explicitly find the wave functions in the 
coordinate representation. In Section III, properties of the squeezed states are analyzed. In the Appendix, we suggest 
a method which allows one to construct symmetry operators for an arbitrary gauge of the vector potential and apply 
the method to the case of a homogeneous magnetic field. 



I. GENERAL PROPERTIES 



H=±-[p--A) , (2) 



The Hamiltonian of a particle with mass m and electric charge e moving in the x — y plane in a magnetic field reads 

1 

2m 

here the vector potential A(A x ,A y ) corresponds to a homogeneous magnetic field perpendicular to the plane, 
(rot A) z = B. The choice of signs in some of the below formulae depends on the sign of eB\ for definiteness, 
we assume eB > 0. In terms of the ladder operators 

d =~f=- 1 (Vx+Wy) , c< = -jJ- — (v x -iv y ) (3) 

V ZlLOc V 21UJ C 

where v x ,y are the non-commuting components of the velocity operator 

h 2 
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p--A, [v x ,v y ] = i— ^ i ( 4 ) 
c mH A 



(the cyclotron frequency lu c = and the magnetic length I = y/hc/\eB\), the Hamiltonian Eq.(||) can be conve- 
niently written as 

H = hu> c (^c+^J , [c,ct]=l. (5) 

In the presence of a homogeneous magnetic field, the translations in the x-y plane and rotations around the z— axis 
remain symmetry elements. The reflection, a v , in a plane passing through the z— axis ( yOz, for definiteness), reverses 
the magnetic field and is not a symmetry transformation. However, the product o~t = T ■ a v of time-reversal T and 
a v , both reversing the field, is a valid symmetry. 

It is well-known that the Hamiltonian Eq.(^) may not commute with the operators associated with the physical 
symmetries because the vector field A(r) has a lower symmetry than the corresponding magnetic field. If this is the 
case, the coordinate transformation should be accompanied by a certain gauge transformation which compensates the 
change in A(r) ||. A procedure of constructing the transformation, which is valid for arbitrary gauge of the vector 
potential, is presented in the Appendix. 

As shown in the Appendix, the operators of finite translations, r — * r + a, are built of the generator of translations 

P = p--A+-B X r. (6) 

c c 

The operator P commutes with the Hamiltonian Eq. (|2|) as manifestation of the translation invariance preserved in a 
homogeneous magnetic field. The components of P obey the commutation relation 

[P x ,P y ]=-£hB . (7) 
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Equation (||) is valid in an arbitrary gauge of the vector potential A. In case of the symmetric gauge, A = ^(B X r), 
Eq.(||) gives the expression first found by Zak ||. 

Presenting P in the form in Eq. ([!]), one indeed recognizes in R the center of the Larmor circle (the guiding center), 



~ mc „ 

R = r + — 2 vXB, 

an integral of motion known from classical mechanics. The following commutation relations 



(8) 



(A), . (*), 



= ihSi, 



- - I 2 
[X,Y] = L 



(9) 
(10) 



can be readily derived from Eqs. ([?]), and (^J). 



II. SQUEEZED STATES 

From Eq.([To|), the X- and Y"-coordinates of the Larmor circle center are incompatible quantum variables. Given 
the commutator Eq. (|To|) , their variances obey the standard (see e.g. uncertainty relation: 

(AX) 2 (AY) 2 > , 

where the variance of a variable A is defined as (A A) 2 = ( I A A) ) with A A here and below standing for A A = 



A — {^Aj. If the uncertainty relation is satisfied with the equality sign, it is said that the system is in a minimum 

uncertainty state or, in other words, in a coherent or, more generally, in a squeezed state. 

We construct the stationary squeezed state, \R, N) as a simultaneous eigenfunction of the Hamiltonian Eq.(^J) and 
the operator 

X* =i"cos$ + Ysin$ , (11) 

where $ is a complex parameter $ = <&i +i$2- Under a ^-rotation around the z-axis, $ transforms to $ — > $' = — 4>, 
and X$ can be also represented as 

1$ = X' cosh $ 2 +if'sinh$ 2 (12) 

where X' = Xcos$i + Y sin$i and Y' = —XsinQi + Y cos$i are the Cartesian components of the guiding center 
R in the principal axes where $' x = 0. 

The state \R, N) is found from the following system of equations 

H\R, N) = hw c (N + ^)\R,N) , (13) 

X 9 \R,N) = (Xoos<f> + YBm$)\R,N) . (14) 

The quantum numbers of a state are the Landau level number N, and the expectation value of the guiding center 
position, R{X 1 Y)] two real parameters X and Y specify the complex eigenvalue X$. 

The operator X$ is not Hcrmitian and the eigenvalue problem may or may not have solutions among physically 
admissible normalizable functions, (R, N\R, N) = 1. To find necessary conditions for the existence of physical 
solutions, we note that Eq.([l4|) leads to (R, N\AX$AX$\R, N) = 0, or, using the representation in Eq.(|l2|) and 
Eq.©, 

(AX 1 ) 2 + tanh 2 ^(AY"') 2 = —I 2 tanh$ 2 • (15) 

Observing that the l.h.s. is positive definite, we conclude that Eq.([l4j) has normalizable solutions only if $2 < (or, 
more generally, eB$>2 < ). 

The real and imaginary parts of the relation (R, N\ (^AX^j \R, N) — 0, give 
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(AY') tanh" $ 2 = (AX ) (16) 
(AX' AY' + AY'AX'^ = (17) 

(AX') 2 (AY') 2 = l - . (18) 



The last relation follows from the first two combined with Eq.([L5|). Also, 

(AA') 2 = ^|tanh$ 2 | , (Ar') 2 = ^|coth$ 2 | . (19) 

According to Eq.(|lg|), the eigenfunctions of X<& indeed belong to the class of minimum uncertainty states. From 
Eq. (|l7|) , one sees that the physical meaning of <&i is to show the orientation of the principal axes, along which the 
quantum fluctuations of the guiding center position are independent. It follows from Eq.(|l^) that <&2 controls the 
relative uncertainty of the projection of the guiding center onto the principal axes. 

To find the wave functions of the squeezed states, we first consider the states from the Landau level N = 0, \R, 0) 
and solve the following system of equations 

c\R,0) =0 , (20) 
X^\R,Q) = X$\R,0) . (21) 

(Eq.(p0[) is equivalent to Eq.(|T3]) for the ground state N = 0.) To find the explicit form of the wave functions, we 
choose the symmetric gauge 

A = ±Bxr 7 (22) 

where different directions are treated on equal footing. 
In the notations 

= (x - X) cos $ + (y - Y ) sin $ , 
yq, = -(x - X) sin $ + (y - Y) cos $ , (23) 

(y^ = 5$+i), the operators in Eq.(^o|) and Eq.(^ll) take the form 

\/2 d . . 8 , 1 . . .. x . 1 



i—r<i " l *c = — h I 



d 

A$ = x* -H 2 ^- + X* , (25) 
dy<s> 

where Y$ = —X sin <f> + Y cos 

In the coordinate representation, Eqs.( [20"|]2l| ) become a system of linear differential equations of the first order for 
^(r\0,R) = (r\R, 0). The normalized solution reads 



*(r | 0, R) — C$ exp f — (^1 + ix<s>y<s> + ix^Y® - iy$A$)^ , 



(26) 



where C$ is the normalization constant 



\C*\ 2 = ^jaV^ -exp 2z($* - I) . (27) 

As expected, the function in Eq.(p6|) is normalizable and Eq.([2^) is meaningful, only if 7m<I> < 0. 

The normalized states for the N-th Landau level, eigenfunctions of the Hamiltonian in Eq.(|5|), can be now found 
from 

*(r\N,R) = -j=(&) rr *(r\0,R). (28) 



The calculation can be easily done with the help of the following identity 
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ct = e- i *e A (-^=-^-) e" A , 
\iV2 dx ® J 

A = — j(x% - ix<s,y$ - ix^Y® + iy$X$) . (29) 

After some algebra we obtain, 

9(r\N, R) = C *^7^=- ex P (*2^ S ■ (-R x r) — ^ai* (x# + ij?*)) ^ (^fj (30) 

where ffjv(0 is the Hermite polynomial, i/jv(£) = (— l) jv e^ 2 d — . In the coordinate representation, this expression 
gives the wave function of the squeezed state \R, N) centered at R and belonging to the iV-th Landau level (in the 
symmetric gauge Eq.(|2])); the "rotation angle" $ in Eq.([2^) is a complex parameter, Im& < 0. 



III. PROPERTIES 

The basic features of squeezed states can be seen in Figs. [I] and ||, where the density and the current are plotted 
for a typical state: N = 2, 5> = i$2, | tanh$2| = 0.1. The squeezed state is localized in the sense that it has a finite 
extension in both x— and y— directions. Qualitatively, the squeezed state is a superposition of classical Larmor orbits 

of radius = J ; the centers of which are positioned in the vicinity of R with a typical deviation AX' and 

AY"'. From Eq. (|l9|) , AX'/ AY' = | tanh<!>2| < 1 so that the state is elongated in the direction of the principal j/'-axis. 
When $2 tends to zero, the elongation increases and the squeezed state asymptotically transforms into a "strip" (of 
length ~ Z/|$ 2 |). 

The wave function of a squeezed state from the N-th Landau level has N isolated zeroes in the x — y plane. 
The zeroes are at the points on the line Im x$ = 0, where the Hermite polynomial H^i^f-) has its N roots. In 
the limit $2 — > —00, the zeroes gather together at the point r — R. This limit, where AX — AY = |, gives the 
stationary coherent state introduced by Malkin and Man'ko fij] , the angular momentum eigenstate with the eigenvalue 
L z = —HN. In the coherent states, the probability density is rotationally invariant. 

In quantum optics, squeezed states can be presented as the result of the action of the "squeezing operator" on the 
coherent state. Similar to squeezing of the cylindrically symmetric coherent states (7m$ — > —00) is achieved by 
applying 

S - exp ^r{x'Y' + Y'X') J , (31) 

where tanh r = e 2 * 2 . In optics there are certain non-linear processes with the evolution operator in the form of S [[l6| . 
A coherent state then evolves into the squeezed state. For the case of a particle in a magnetic field, the analogous 
problem of preparation of a squeezed state has been considered in [|l3| . 

As discussed in the Appendix, the product of the mirror and time reversal transformations is a valid symmetry 
element. As a consequence, the distribution of the density and current are mirror symmetric (relative to the principal 
axes) as also apparent in Figs.(Q-^). 

Within a given Landau level, the squeezed states, eigenfunctions of a non-Hermitian operator A"$, are non- 
orthogonal. For the states defined in Eq.([30|) with the real positive normalization constant C$ Eq. (p7|) , the overlap 
integral reads 

\ 1/4 / \ 1/4 

1 - exp (2i($* - $)) ) 1 1 - exp (2i($'* - $')) 
(R; $, N I R'; N') = 5 NN , ± ; '- ^ — '— x 



(l - exp (2i($* - $')) 



/ \ 1/2 

(Re I 1 — exp (2i($* — $')) j > 0). The overlap of the states differing in the Larmor center position R or the 
parameter $ does not depend on the Landau level number N as it follows from the Eq.d28) and the commutation 
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relation in Eq .(p[) - Therefore, the overlap integral in Eq.(|32|) can be calculated using the Gaussian wave functions for 
N = in Eq.(H). 

Repeating the derivation known in the theory of coherent states of a harmonic oscillator (see e.g. ||), one can show 
that the set of squeezed states is complete i. e. the closure relation, 

°° r rP-Ti 

i = E Y%\R,N)(R,m, (33) 

is valid. As in the harmonic oscillator case, the states \R, N) with continuously varying R form an overcomplete set 
within the N-th Landau level. Repeating Perelomov's arguments fL7|] , one can show that the subset of states with 
i?'s on the sites of a periodic lattice is overcomplete if the lattice is too dense (the unit cell area so < 2irl 2 ) and is 
not complete for a too dilute lattice ( sq > 2irl 2 ). When sq — 2tt1 2 (i.e. the flux through the unit cell equals to the 
flux quantum — ), the system of the functions is complete and it remains complete even if a single state is removed; 
it becomes incomplete, however, if any two states are removed. 

In conclusion, stationary squeezed states of a charged particle in a homogeneous magnetic field have been constructed 
and analyzed. The distinctive feature of the squeezed states is the minimal quantum uncertainty of the position of the 
Larmor circle center. The family of the squeezed states is characterized by the squeezing parameter Im $ variation of 
which allows one to transform gradually strip-like states (infinitely extended in the direction controlled by Re $) to 
eigenfunctions of the angular momentum with rotationally invariant density distribution. The squeezed states have a 
rather simple structure: As it follows from Eqs.([50] and |23|), a general squeezed state can be obtained from a Landau 
"strip" by a complex angle rotation of the coordinates. The simplicity of the construction gives the hope that the 
squeezed states may turn out to be useful. 

We are grateful to J0rgen Rammer for reading the manuscript and valuable remarks. We also thank V.V. Dodonov 
for drawing our attention to the references |12|-|l4|l . This work was supported by the Swedish Natural Science Research 
Council. 



APPENDIX A: SYMMETRIES OPERATORS 



The physical symmetry of a system in an external magnetic field is controlled by the symmetry of the vector field 
B(r) (among other factors). However, the Hamiltonian, H — Ha contains the vector potential A(r) rather than B. 
The vector field A(r) has a lower symmetry and, moreover, the spatial symmetry of A(r) is gauge dependent. For this 
reason, the Hamiltonian often does not commute with the operator corresponding to a physical symmetry element. A 
homogeneous magnetic field gives a simple example: the translation invariance is preserved but the vector potential is 
always r-dependent. It is well-known |3|] that the symmetry operator should include a certain gauge transformation 
which compensates the change of the A-field generated by the symmetry transformation. 

To build the modified operators on a regular basis, we suggest the following procedure. For any given symmetry 
element O, it is always possible to find field A(°'(r), rot A^ = B, which is invariant relative to O: OA^ — A^°\ 
Whatever gauge is chosen for A in Eq.(||), the gauge transformation, Gq 1 HaGo, specified by 

r 

Go = e^ 0) , x <°>(r) = ±Jdr(A- A<°>) , (Al) 



changes the vector potential entering the Hamiltonian from A to A — VrX^ = A^°\ i.e. Gq^HaGq = H A (o) . 
By construction, O does not change the vector field A^°\r) and, therefore, the transformed operator Gq 1 HaGo 
commutes with O, i.e. 

(g^HaGo) 6 = 6 (g^HaGo) or H A [g 6G 1 ) = (g 6G^) H a ■ 

Therefore, it is the operator 

6 a = g 6g 1 , (A2) 

which commutes with the Hamiltonian Ha and represents the symmetry clement O. As such, the operator Oa depends 
on the gauge chosen for the vector potential A in the Hamiltonian, but its matrix elements are gauge invariant if the 
sandwiching functions are gauge transformed in the usual manner. 



G 



Equivalently, the symmetry operator in Eq.(|A2|) can be written as 

a A = e *(* (0, M-* <0) «M 6 (A3) 

= a e i(x (o) (6-M-x""w), (A4) 

Below, we analyze only the case of a homogeneous magnetic field but the method is generally applicable. 

In the presence of a homogeneous magnetic field, the translations in the x-y plane and rotations around the z— axis 
remain (continuous) symmetry elements. As the reflection, o v , in a plane passing through the z— axis ( yOz, for 
definiteness) , reverses the magnetic field, o v is not a symmetry element. We note here that the product of time- 
reversal T and o v , both reversing the field, is a symmetry element. We denote the product by ot = T ■ o v . 

o 

The transformation O changes the wave function as, ip — > Oip, where O denotes the operator associated with O. 
For translations and rotations, the form of the operator O is obvious; in the case of ot, 

a T ip(x,y) = ip*(-x,y) , (A5) 

and <jt is an anti-linear anti-unitary operator. From 0~ 1 HaO = Hqa, the transformed vector field OA can be 
found. Again, the result is obvious in case of rotations and translations. Under ot, which is the mirror reflection 
of the polar vector field A(r) in combination with time reversal (A — ► —A), the vector potential transforms as: 
o T A x (x,y) = A x (-x,y), o T A y (x,y) = -A y (-x,y). 

First, we consider the operator of a finite translation O = T a , T a tp{x, y) = ip(x + a x , y + a y ), a being the translation 
vector. The vector potential, A^ = A^ a \ invariant relative to the translation along a, A(°'(r) = A^(r + a), can 
be taken as, 

A {a) = n((B X r)-n) , a = \a\n , (A6) 



("Landau gauge"). Using Eq . (|A3|) , the operator of magnetic translation reads 



iJjj j dr'( A(r')-A (<1 '(r') 

f A a) = e C V ' e * p - a , (A7) 

p = jV being the canonical momentum. As a consequence of the physical translational invariance, this operator 
commutes with the Hamiltonian Eq.(||) for arbitrary gauge of the vector potential A(r). 

Up to terms linear in a, f { A a) = 1 + |a- (p - f A+ f A^) or, using Eq.@, f A a) = 1 + |a - (p - |A + f-B X r). 
From here, one reads off the expression for the generator of translations P and derives Eq.(^). 

This derivation of Eq. (^) links P to the translational symmetry. The same expression for P can be derived in a 
more intuitive manner: First, one considers the operator of the guiding center R, the expression for which in Eq.(|^) 
can be guessed from the correspondence principle. Since R is a classical integral of motion, R must commute with 
the Hamiltonian. Now, one defines P by Eq.(Q) and comes immediately to Eq.(^). 

Next we consider rotations i.e. O — R. The "bare" operator of a rotation around the z-axis is R = e^ 7 "*^ 



being the angle of the rotation. The R— invariant vector potential is A^ = x r. Applying Eq.(A2), the symmetry 
operator reads: 

R A {(t>) = e C o e i(rxp),<f> e ° o . (A8) 

In the limit (f> — > 0, Ra(<I>) ~ 1 + i£>z<i>, where the generator of rotations C z is 

C z = e (rxp) z e 

Simplifying and using the expression for A^ , one gets 

£ z = (rx(p-~A)) +-?-(rx(Bxr)) z . (A9) 
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For the case of a homogeneous magnetic field, this operator commutes with the Hamiltonian for arbitrary gauge of 
the vector potential and reduces to the usual angular momentum L z = rxp when A, B — > 0. Also, C z — L z in the 
symmetric gauge A — \E> X r. 

The operator C z can be written in the following identically equivalent forms 

- 1 eB -2 mc « , . s 
= (A10) 

- 1 c ~ 2 mc - , . , 

£ ' = 5S p -S ff! (A11) 

iJ, and the (two-dimensional) vectors P and -R are defined in Eqs.(||), (|J) and (0), respectively. One sees that the 
integral of motion C z is, actually, a function of the other conserving quantities P (or R) and the energy H. 

If the axis of the a rotation is shifted from the origin to the point Rq 7 the generator of the rotations denoted as 
£r ,z reads 

C Ro , z = \^(h-R«) 2 -^H. (A12) 

Note that the vector potential enters C z Eq. (|A9|) and the generator of translations P Eq.(^) only in the gauge 
covariant combination p — |j4 so that their matrix elements are gauge invariant. 

Finally, we consider the combined mirror and time-reversal transformation in Eq.( ]A5[ ): O = <jt- One can check 
that &t does not change the Hamiltonian in Eq.(||) if the symmetric gauge Eq.(p2|) is chosen. Ultimately, this is the 
reason for the mirror symmetry in the distribution of the density and current seen in Figs([l]- Q). 
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FIG. 1. The probability density for the state N = 2 ; tanh<E> 2 | = 0.1, $i = located at the origin R = 0. 
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FIG. 2. The current density for the state iV = 2 ; tanh 3> 2 | = 0.1, <J>i = located at the origin Ji = 0. 
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FIG. 4. The probability density for the state N = 10 ; | tanh $ 2 1 = 0.1, $1 = located at the origin R = 0. 
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